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ABSTRACT 

Context. The nearby universe remains the best laboratory to understand the physical properties of galaxies and is a reference for any 
comparison with high redshift observations. The all sky (or very large) surveys that have been performed from the ultraviolet (UV) to 
the far-infrared (far-IR) provide us with large datasets of very large wavelength coverage to perform a reference study. 
Aims. We investigate the dust attenuation characteristics, as well as the star formation rate (SFR) calibrations of a sample of nearby 
galaxies observed over 13 bands from 0.15 to 160 lum. 

Methods. A sample of 363 galaxies is built from the AKARI /FIS all sky survey cross-correlated with SDSS and GALEX surveys. 
Broad band spectral energy distributions are fitted with the CIGALE code optimized to analyse variations in dust attenuation curves 
and SFR measurements and based on an energetic budget between the stellar and dust emission. 

Results. Our galaxy sample is primarily selected in far-IR and mostly constituted of massive, actively star-forming galaxies. There 
is some evidence for a dust attenuation law that is slightly steeper than that used for starburst galaxies but we are unable to constrain 
the presence or not of a bump at 220 nm. We confirm that a time-dependent dust attenuation is necessary to perform the best fits. 
Various calibrations of the dust attenuation in the UV as a function of UV-optical colours are discussed. A calibration of the current 
SFR combining UV and total IR emissions is proposed with an accurate estimate of dust heating by old stars. For the whole sample, 
17% of the total dust luminosity is unrelated to the recent star formation. 

Key words, galaxies: star formation-(ISM:) dust, extinction-infrared: galaxies-ultraviolet: galaxies 



1. Introduction 

The broad-band spectral energy distribution (SED) of a galaxy 
represents the combination of the emission from both stars of 
all ages and interstellar dust that interact in a complex way by 
means of the absorption and scatteiing of the stellar light by dust 
grains. The stars emit from the UV to the near-lR whereas the 
mid and far-IR emission comes from interstellar dust heated by 
the stellar emission. The broad-band SEDs can be reconstructed 
from a library of stellar tracks by assuming a star formation his- 
tory, an initial mass function, and a scenario for dust attenua- 
tion to reproduce the observed UV-optical distribution and the 
re-emitted emission in mid and far-IR. By comparing data with 
models, one can attempt to derive some physical parameters re- 
lated to the star formation history and dust attenuation in a ho- 
mogeneous way and simultaneously for all galaxies of a given 
sample. Practically the current star formation rate (SFR), the 
stellar mass, and the amount of dust attenuation are commonly 
derived from broad-band photometry (e.g. ,Salim et al. 2007 da 



ICunha et al.| [20081 [Walcher et all |2008l [20T0l and references 
therein). Without spectral information and given the high degree 
of degeneracy in the SEDs, details about the star formation his- 
tory such as stellar population ages, or specific aspects of dust 
attenuation such as its wavelength dependence, are more diffi- 
cuk to constrain (e.g. jWalcher et al.|[2&T0l|Noll et"aLl|2009b| l . 
With the availability of mid and far-IR data for large samples of 
galaxies, codes that combine stellar and dust emission to anal- 
yse SEDs are particularly useful. The standard approach consists 
of solving the radiation transfer in model galaxies to build self- 
consistent SEDs from the UV to the far-IR (e.g. jPopescu et al7 



2000; Sil va et al.| [T9981 jTuffs stU] [2004). These sophisticated 
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models require complex calculations and are not directly appli- 
cable to large samples of galaxies. Libraries of templates can be 
produced ( [Iglesias-Paramo et"ar 2007 [ Popescu et al. 2010 1, as 
well as some recipes ( [Tuffs et al. 20041, but the various free pa- 
rameters and physical assumptions are difficult to constrain with 
only broad-band, integrated fluxes. Another tool applicable to 
large samples of galaxies consists of models based on a simple 
energetic budget where the global dust emission corresponds to 



2 



V. Buat et al.: AKARI-SDSS-GALEX galaxies 



Table 1. Sample of galaxies, data with an asterix (*) are not used 
to fit the SEDs with CIGALE 



Survey 



Band 



Wavelength 



Nb of sources 



GALEX FUV 153 nm 363 

GALEX NUV 231 nm 363 

SDSS u,g,r,i,z 355, 469, 617, 748, 893 nm 363 

AKARI/FIS' N60 65 yum 129 

AKARI/FIS' WIDE-S 90 yum 357 

IRAS Band-3 60 359 

IRAS Band-4 100 yum 354 

AKARl/FlS WIDE-L 140 //m 363 

AKARI/FIS N160 160 //m 111 



the difference between the emitted and observed stellar light (Idal 
[Cunha et al.][2008j|Noireta l. 2009b). Codes based on synthetic 
stellar population modeling and simple emission properties for 
the dust are well suited to analysing large datasets at the expense 
of an over simplification of the physical processes at work in 
galaxies. 

The natural laboratory for this SED analysis is the nearby uni- 
verse for which plenty of high-quality photometric data are 
available. 

Several studies have been based on the Spitzer Nearby 
Galaxies Survey (SINGS) ( |Kennicutt| |2003| l, which provided 
broad-band spectral energy distributions from 0.15 to 160 yum 
for 75 galaxies pale et al.| [2007) 1. |da Cunha etal] ( |2008| and 
Noll et al. ( 2009b I illustrate the capabilities of their fitting code 



with this sample. 

In the present work, we investigate the properties of nearby 
galaxies by selecting a sample of objects that is as large as pos- 
sible, observed during the AKARI, SDSS, and GALEX surveys 
that correspond to 13 photometric bands from 0.15 to 160 yum. 
Particular care was taken to measure total flux densities for the 
extended sources. Our aim is to perform SED fitting to deduce 
physical properties in a homogeneous and consistent way, focus- 
ing on dust attenuation and SFR measurements. The originality 
of our study is to extend the wavelength coverage beyond 100 
jum (the longest wavelength observed by IRAS) so that we con- 
sider only galaxies with a reliable flux at 140 //m: our selection 
is close to a selection at 140 fim and consists of 363 unresolved 
objects with a very large wavelength coverage and secure mea- 
surements, allowing a statistical analysis. A selection at 140 yum, 
beyond the peak of dust emission for most galaxies is particu- 
larly valuable for further comparison with distant galaxies ob- 
served in sub-mm by Herschel. It is also the first application of 
our code CIGALE to a large sample of nearby galaxies with such 
a wavelength coverage. 

In section 3, we describe the SED fitting method, the reliabil- 
ity of its results derived by the analysing a mock catalogue and 
the specific results obtained for our galaxy sample. Section 4 
is devoted to the analysis of dust attenuation characteristics and 
section 5 to SFR calibrations. The conclusion is presented in sec- 
tion 6. All magnitudes are given in AB system. We assume that 
Q,„ = 0.3, Qa = 0.7, and Hq = 70 km s"' Mpc"'. 

2. Data 

Our aim is to build a galaxy sample with high quality fluxes 
from the UV to the far-IR. The present sample is built from the 
primary sample of |Takeuchi et al.|p010| l. It consists of galax- 
ies observed as part of the AKARIIFIS all sky survey (Kawada 
|et al.| |2007|l and restri cted to the area covered by SDSS/DR7 
( [Abazajian et al.| |2009| ) and the GALEX All Imaging Survey 
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Fig. 1. AKARI fluxes at 140 yim plotted against IRAS fluxes at 
100 [im (upper panel) and AKARI fluxes at 90 yum (lower panel). 
The sample is divided into three sub-samples according to the 
diameter d of the galaxies: d < 1' (dots), 1' < d < 1.5' (small 
circles), and d > 1.5' (large circles). The data are compared to 
model predictions from Dale & Helou ( 2002| l for three different 
values of the a parameter, a — 1.5,2, and 2.5 (dotted, solid, and 
dashed lines respectively). 



( Martin et al.| |2005| l for an additional cross-match. [Takeuchi] 
|et al. (2010 1 cross-matched the first primary catalogue of the 
AKARll^l^ All Sky S urvey (j6-l version) with the //M5/PSCz 
( Saunders et al. 2000[ ) to select 776 sources confirmed as galax- 
ies, with a known redshift and observed with GALEX. For 607 
of these 776 sources, SDSS images are suitable for perform- 
ing accurate photometry (no bright star superposed). We refer 
to |Takeuchi et al.|pOIO] ) for details of the sample selection and 
photometric measurements. Briefly, specific aperture photome- 
try was performed for each source in the GALEX and SDSS im- 
ages at the AKARI positions to avoid shredding. Total flux den- 
sities were measured in elliptical apertures on the GALEX im- 
ages using the method described in Iglesias-Paramo et al. ( 2007[ l. 
When sources were shredded in the SDSS images during the 
deblending step of the pipeline (most of the cases), we used 
the photometry of the parent object measured before deblend- 
ing within the Petrosian radius. The errors due to the photomet- 
ric measurements at UV and optical wavelengths were measured 
and are of the order of ~ 5%. 

We collect the fluxes from the A KARI/FIS All-Sky Survey 
Bright Source Catalogue Version 1.0 (Yamamura et al. 2010 1, 
whose errors are smaller than those of the /3-1 version used in 
[Takeuchi et al.| ( [20T0| ). We restrict the sample to galaxies whose 
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Parameters 


Symbol 


Range 


Star formation history 


metallicities (solar metallicity) 


Z 


0.02 


T of old stellar population models in Gyr 


Tl 


1; 3.0; 5.0; 10.0 


ages of old stellar population models in Gyr 


h 


13 


ages of young stellar population models in Gyr 


h 


0.025; 0.05; 0.1 ;0.3 ;0.5 ;1.0 


fraction of young stellar population 


fvSP 


0.001; 0.01; 0.1; 0.999 


IMF 


K 


Kroupa 


Dust attenuation 


Slope correction of the Calzetti law 


S 


-0.3; -0.2; -0.1; 0.0; 0.1; 0.2 


V-band attenuation for the young stellar population 




0.15; 0.30; 0.45; 0.60; 0.75; 0.90; 1.05; 1.20; 1.35; 






1.5; 1.65; 1.8; 1.95; 2.1 


Reduction of Av basic for old SP model 


/all 


0.0; 0.50; 1.0 


IRSED 


IR power-law slope 


a 


1.0; 1.5; 1.75; 2.0; 2.25; 2.5; 4.0 



Table 2. List of the input parameters of the code CIGALE and their selected range. 



140 yum flux is of high quality (FQUAL140=3), 363 sources 
fulfiling this condition. We find that most of selected sources 
(357 sources) also have high quality fluxes at 90 fj.m but the 
fraction of these sources with reliable fluxes at 65 and 160 
yum is about ~ 1/3. The fluxes uncertainties are estimated to 
be ~ 20% as quoted in the AKARI/FIS AU-Sky Survey Bright 
Source Catalogue. All the 363 sources were detected by GALEX 
and in all SDSS bands. 

In Table 1, we gather the number of sources with high qual- 
ity fluxes at each wavelength of interest, the sample used for 
our SED analysis consisting of the 363 sources reliably detected 
at 140 fj.m and in all the SDSS bands. Thirteen bands are con- 
sidered, most galaxies being detected in 11 of them. All the 
fluxes are corrected for Galactic extinction using the Schlegel| 
|et al.| ( |1998l l dust maps and thc jCardelU et aI.| ( [TM9| extinction 
curve. The redshift distribution of the sources is similar to that of 



the original sample of |Takeuchi et al. (20101 with a mean value 
<z>= 0.016 

Since the detection rate at 65 /jm is low (129/363 sources), we 
prefer to consider the IRAS fluxes at 60 fj.m available for 359 
sources (we use the 60 fim co-added fluxes from the PSCz and 
consider only sources with an IRAS quality flag equal to 3). 
AKARI fluxes at 65 and 90 yum are systematically lower than 



IRAS fluxes at 60 and 100 //m, respectively (Takeuchi et al. 



2010] ). This discrepancy is likely to be due to the higher spatia 
resolution of AKARI data (pixel size 27 arcsec, FWHM of the 
PSF equal to 39 arcsec ( [Kawada et al.|[2007l l): the point source 
extraction is no longer valid for slightly extended nearby sources 
(Yuan et al. in preparation). Therefore, at 100 fim we also con- 
sider IRAS co-added fluxes for the 354 sources safely detected 
by IRAS at 100 fim (IRAS quality flag equal to 3), in section 3 
we discuss the difference in the estimate of the total IR luminos- 
ity LiR when using either only AKARI data or a combination of 
AKARI and IRAS data. Optical diameters measured at the isopho- 
tal level B= 25 mag arcsec"^ are compiled from the HyperLeda 
database Q 

Any limitation of the spatial resolution and point source extrac- 
tion are expected to be less important at 140 and 160 fim than at 
65 and 90 /vm, with a pixel size of 44 arcsec and a PSF FWHM of 



58 and 61 arcsec, respectively, at 140 and 160 jim ( Kawada et al. 
2007| l. In Fig[T[ the 140 //m fluxes are compared to the fluxes at 



90 and 100 fim together with the predictions of the |Dale & HeIou| 
(2002 ) models representative of galaxies in the nearby universe 
(corresponding to an a parameter of the Dale & Helou ( 2002| l 



models between 1.5 and 2.5 as described in section 3) . Models 
and IRAS data at 100 yum are found to be consistent without any 
obvious trend with the galaxy diameters whereas AKARI fluxes 
at 90 fim appear systematically underestimated, the discrepancy 
increasing with the optical diameter, confirming the analysis re- 



sults of Takeuchi et al. ( 2010 1 and Yuan et al. (in preparation). 



3. SED fitting 

3.1. CIGALE code 

An efficient way to derive physical parameters of star forma- 
tion and dust attenuation homogeneously is to fit the observed 
SED with models from a stellar population synthesis code. We 
use the code CIGALE (Code Investigating GALaxy Emission) 
which derives physical information about galaxies by fitting 
their UV-to-far-IR SED ( |NoU et al.[ |2009bt [Giovannoli eFaL] 
|2010 l. A Bayesian analysis is used to derive galaxy properties 
similar to that developed by Kauffmann et al. (2003a). CIGALE 
combines a UV-optical stellar SED and a dust, IR-emitting com- 
ponent. First, models are built, then each model is quantitatively 
compared to the observed SEDs as we account for the uncer- 
tainties in the observed fluxes. The probability function of each 
parameter is calculated and the estimated value of the parameter 
and its error coiTespond to the mean and standard deviation of 
this distribution. Models are generated with a stellar population 
synthesis code, assuming a particular star formation history and 
dust attenuation scenario. The energetic balance between dust- 
enshrouded stellar emission and re-emission in the IR is care- 
fully conserved by combining the UV-optical and IR SEDs. We 
refer the reader to [Noll et al. (2009b I for a detailed description. 
Briefly, the IR SEDs are built from flie |Dale & Heloul ( [2002 ) tem- 
plates. CIGALE allows for additional dust emission from a dust- 
enshrouded AGN but does not include the unobscured emission 
of an active nucleus. Therefore, we discard five galaxies clas- 
sified as Seyfert 1 by Veron-Cetty & Veron ( 2006[ l. We then 
performed several tests after adding two AGN templates (corre- 
sponding to a PAH-free emission with respectively Ay = 32 and 
64 mag of extinction from the |Siebenmorgen & Krii gel (2007| 
library) to the SED fitting (Noll et al.||2009b| l and found that the 
AGN contribution to the total IR emission never exceeds 12% 
(without any difference for galaxies classified as active or not) 
but is very badly constrained in the absence of mid-IR data (a 
null contribution of the AGN component is never excluded by 
the Bayesian analysis). 



' http://leda.univl-lyonl.fr/ 



^ http://www.oamp.fr/cigale 
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Fig. 2. Scenario of star formation history adopted in this work. 
The old stellar population is produced by an exponentially de- 
creasing star formation rate. By following a steep-to-moderate 
e-folding rate ti, the young stellar population is created in t2 
years, at a constant, adjustable rate. 
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Fig. 3. FUV-NUV versus NUV-r colours for observed (black 
point) and modeled (grey point) data. 



In the following, we describe the parameters that are crucial 
to this study. The input parameters values used in this work are 
listed in Table|2]and presented in this section. They were chosen 
after several trials and represent a compromise between a good 
estimate of the parameters of interest and the duration of the 
SED fitting process. An over sampling of ill-constrained param- 
eters may also reduce the accuracy of well-constrained parame- 
ters by increasing the number of models of moderate probability 
and consequently the dispersion in the parameter estimates. This 
effect was found to be moderate but real, thus there is no need 
to perform a rigorous sampling of all the parameters. The range 
of values was chosen to be large for the first tests, then reduced 
by excluding values never chosen during the minimization. 
The sampling of the parameters was then optimized in terms of 
computer time and parameter estimation. 

3.1 .1 . Stellar populations and star formation histories 



Fig. 4. FUV-NUV versus NUV-r colours for observed (black 
point) and mock (empty square) data. 



asymptotic giant branch stars. The chosen initial mass function 
is that of Kroupa^(2001 ). The metallicity is taken to be solar. 

The star formation history implemented in CIGALE is the 
combination of two stellar components resembling an old and 
a young stellar populations. The two populations roughly rep- 
resent a burst of star formation in addition to a more passively 
evolving stellar component. We verify that including these two 
populations improves the results of the SED fitting relative to 
the case where only a single exponentially decreasing star for- 
mation rate is assumed. The reality is certainly more complex 
but adding more components may lead to degenerate solutions. 
The old stellar component is modeled with an exponentially de- 
creasing SFR (with various values of the e-folding rate ti) that 
started 13 Gyr ago (fi - 13 Gyr), the young stellar component 
consisting of a burst of constant star formation starting later ( f2 
Gyr ago) whose amplitude is adjustable (Fig.[2]i. We adopt four 
values of ti (1, 3, 5, or 10 Gyr) and is taken in the range 0.025 
to 1 Gyr. The two stellar components are linked by their mass 
fraction, /ysp, which corresponds to the ratio of the mass locked 
in the young stellar population mass to the total stellar mass. The 
parameter /ySp varies in the range 0-1, with a logarithmic scale 
(linear or logarithmic variations are allowed for the input param- 
eters). CIGALE estimates the SFR, defined as (l-/ySp)*SFRi + 
/ySp*SFR2, where SFRi and SFR2 correspond to the star forma- 
tion rates of the old and young stellar populations, respectively. 
In the following, the SFR refers to the formula calculated for the 
current time. 

In its newest version, CIGALE also estimates the dust-free 
D4000 break measured in the spectrum of the composite stellar 

{T999| 



population. The calibration used is that of Balogh et al 



(see also section 3.4 for a more detailed discussion). 



3.1.2. Dust attenuation and infrared emission 

To model the attenuation by dust, the code uses the attenuation 
law of |Calzetti et al. (2000) as a baseline, and offers the pos- 
sibility of varying the steepness of the attenuation law and/or 
adding a bump centred at 220 nm (see INoU et al. ( 2009b^ for 
a complete description of the dust attenuation prescription). To 
produce different s lopes, the parameter 6 is introduced and the 



We adopt the stellar population synthesis models of Maraston 
( |2005| l, which include a full treatment of the thermally pulsating 



original |Calzetti et al. (2000) law is multiplied by (A/AyY with 



Ay = 550 nm. We choose 6 between -0.3 and 0.2. This range of 
values was chosen after several trade-offs and we exclude larger 
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Fig. 5. Results of Bayesian analysis of the mock catalog. Each plot corresponds to an output parameter considered in this work and 
discussed in the text. The first six output parameters are related to the star formation history (Mjtar, SFR, D4000, ti, f2, and /ysp); 
the next five output parameters characterize the dust attenuation (6, Ay.ysp, /att, ^v, and Apuv) and the two last ones define the IR 
dust emission ( Lir and a). The 'exact' values are plotted on the x-axis, the results applying the SED fitting method to the mock 
data are plotted on the y-axis, and the standard error given by the Bayesian analysis is over plotted as an error bar for each value. 
The Unear Pearson correlation coefficient (r) is indicated on each plot. 



departures from the Calzetti law since these solutions were never 
chosen by the SED fitting process. Several tests are made to 
check whether we are able to constrain the presence of a UV 
bump with the available photometric broadbands, including the 
analysis of mock catalogues, described in section 3.3, and we 
conclude that the bump is unconstrained. Thus, we decide to 
only consider slope deviations provided by the 6 parameter 

The stellar population models must be attenuated before the 
IR emission can be added, since the total IR luminosity is defined 
as the dust-absorbed luminosity of the stars. The code allows us 
to consider that the old stellar population is less affected by dust 
attenuation than the young one as recommended by several stud- 
ies (e.g. [Calzetti et al.] [20001 [Chariot & FalU [2000| [Panuzzo et| 



2007| ). The primary input parameter is the dust attenuation 



in the V band of the young stellar population, Av,ysp- A reduc- 
tion factor /att of the dust attenuation for the old stellar popu- 
lation (relative to the young one) is also introduced as an input 
parameter The global dust attenuation for the total (young and 
old) stellar population is defined as the output parameter Ay. The 
range of values adopted for Ay.ySp and /a,t are listed in Table l2| 
Since we are particularly interested in dust attenuation in the ITV 
range, we also define Apuv, the dust attenuation in the FUV fil- 
ter 

To fit IR observations, CIGALE uses the semi-empirical one- 



parameter models of Dale & Helou ( |2002) . These 64 models 
are parametrized by a, which relates the dust mass to heating 
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Fig. 6. Five examples of best-fit models with different reduced minima chosen to be representative of our sample. The observed 
data are plotted with red crosses and the best-fit model with a solid line. The distribution of the reduced minima for the entire 
sample is plotted in the last panel. 



intensity. The parameter a is directly related to 7?(60, 100) = 
fmnmlfmiim, where /eo^,,,, and fim^im are fluxes in Jy at 60 and 
100 //m, respectively. We choose a in the interval [1 ; 4] to cover 
a large domain of activity. 



3.2. Reliability and distribution of the output parameters 

A first check of the ability of the models to reproduce the data is 
to compare their photometric distributions. The parameters we 
wish to constrain are related to the star formation history and 
dust attenuation characteristics that mainly affect the UV emis- 
sion. Therefore, we choose to compare the distributions of the 
FUV-NUV and NUV-r colours, both being very sensitive to the 
star formation history and dust attenuation. One can see in Fig[3] 
that the range of the modeled colours covers the range of the 
observed ones confirming that the stellar populations and the 
dust attenuation assumed in the models clearly reproduce the 
observed characteristics of the sources . 

The reliability of the parameter estimation is then tested by ap- 
plying the method developed by GiovannoU et al. ( 2010| l, which 
consists of building a specific mock catalogue for each set of 
data to be fitted. The mock catalogue is made of artificial galax- 



ies for which we know the flux in each photometric band of inter- 
est and the exact values of the parameters used by CIGALE. To 
construct this catalogue, we follow three steps. The first step is 
to run the code on the data to choose the best-fit model for each 
object by means ot a simple minimization method. We se- 
lect the 28 galaxies detected at all IR wavelengths: 60, 100, 140, 
and 160 fj.m. In the second step, we add a randomly distributed 
error (which represents typically 10% of the flux) to each flux 
of the best-fit models measured in the photometric bands of the 
dataset: we have a catalogue of artificial galaxies which are ob- 
served in the same filters as the real galaxies and representa- 
tive of them. The representativeness of our mock catalogue is 
checked by comparing the colours (FUV-NUV and NUV-r) of 
the simulated and real galaxies (Fig |4]l. The last step is to run 
the code on these simulated data and compare the exact values 
of the parameters with the values estimated by the code. Fig |5] 
compares the output parameters from the Bayesian analysis for 
the mock galaxies considered in this work with the values esti- 
mated by the code. Values of the linear Pearson correlation coef- 
ficient (r) are also indicated in the plots. The output parameters 
for the real data are all Bayesian estimates. We find very good 
correlations for Mstar, SFR, D4000, Apuv, ^v, ^v.ySP. ^ir. and a 
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Fig. 7. Distribution of the Bayesian estimates of the output pa- 
rameters discussed in this work: M^^.^^, SFR, D4000, Apuv^ 
Av,ySP, Uk, a, 6, and /an. 



with r > 0.80. The mass fraction /ysp is only poorly estimated (r 
= 0.71) and the estimations of ri, t2, 6, and /an are much less sat- 
isfying, with a correlation coefficient equal to 0.15, 0.52, 0.54, 
and 0.45 respectively. According to this analysis which was also 
confirmed by that of the Bayesian error distribution ( |Noll et aL] 



2009b I, we conclude that the reliable parameters are Mjtar, SFR, 
D4000, Apuv, Ay, Av,ysp, ^ir, and a. These parameters are used 
in the following, with some discussion about S and /att- 

The SED fitting and the Bayesian analysis are performed on 
the dataset with the input parameters listed in Table|2] We restrict 
the analysis to objects for which minimum value of the reduced 

is lower than 10: 342 galaxies satisfy this criterium, the re- 
duced distribution for the entire sample is presented in Fig 
|6] (last panel). The discrepant cases are likely to correspond to 
misidentified objects and GALEX data that are inconsistent with 
the other ones. Given the low number of these catastrophic cases 
(~ 6%), we omit them from any further analysis. Some exam- 
ples of best-fit models superposed on the observed fluxes are 
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log(L,„)(IRAS-AKARI)(Le) 



Fig. 8. Total IR emission Lir obtained by running CIGALE with 
IRAS and AKARI data (x-axis) and with AKARI data only (y- 
axis). The solid line correspond to equal values on both axes. 



given in Fig|6] Fig. |7] shows the distributions for the parameters 
of interest estimated with the Bayesian analysis for the 342 ob- 
jects: Mstar, SFR, D4000, Apuv, Ay, Ay.ysp, LiR, a, 6 and /att- 
Galaxies are massive with a median value of Mstar equal to 
10"'^Mo. We find that the SFR is between 0.1 and 100 Moyr"' 
with very few values higher than 30 M^yi^^ and a median value 
of 4 Moyr"'. The galaxies are moderately luminous (< Lir >= 
10.7 + 0.3 Lq) as expected for a flux-limited sample of nearby 
galaxies. The D4000 distribution corresponds to that of actively 
star-forming galaxies in the nearby universe with a median value 
equal to 1.3 (Kauffmann et al. (2003a) and discussion below). 
The parameter Ay.ysp spans the entire range of input values be- 
tween 0.5 and 2.1 mag with a few objects above 1.25 mag. For 
most of the objects. Ay < 1 mag. These output parameters are 
discussed in more detail in the following sections. 



3.3. Dust emission 

3.3.1. Total infrared luminosities 

As discussed in section 2, AKARI measurements at 65 and 90 
IJ.m seem to be underestimated relative to IRAS fluxes at similar 
wavelengths. We can compare the measure of the total IR lu- 
minosities LiR obtained when considering the four AKARI-FIS 
bands or a combination of IRAS (at 60 and 100 /im) and AKARI- 
FIS (at 140 and 160 jum) data. We run CIGALE in both cases 
and the resulting values of Lir are compared in Fig. [8] Both es- 
timates correlate very well, the values obtained with the four 
AKARI bands being systematically lower than those obtained 
by combining AKARI and IRAS data by 25%. In the following, 
we only consider the combination of IRAS and AKARI data. 

3.3.2. Infrared SEDs 



The IR dust emission is fitted using the [Dale & Helou| ( |2002| l 
library (hereafter DH02). These models are parametrized by a 
value of a that is directly related to R(60,100) the ratio of the 
fluxes at 60 to 100 fim observed by IRAS. In Fig. |9] the val- 
ues of a obtained for galaxies observed at 60, 100, 140, and 
160 fim are plotted against the ratio of the IRAS fluxes at 60 
to 100 fim, R(60,100), and compared to the relation between 
a and R(60,100) from the DH02 models. A very good agree- 
ment is found between a values found with CIGALE and the 
values expected from R(60,100). The SED library of DH02 is 
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R(60,100) 



Fig. 9. a parameter estimated with CIGALE as a function of the 
60 to 100 fim flux ratio, R(60,100) for galaxies observed at all 
far-IR wavelengths, up to 160 yum. Vertical error bars represent 
the standard deviation in the determination of a given by the 
Bayesian analysis. The relation between a and R(60,100) for the 
DH02 models is plotted as a sohd hne. 



built to clearly reproduce the far-IR (longward 100 fim) and sub- 
mm (at 850 ;um) emission of local galaxies. Therefore, good 
agreement is expected when data longer than 100 fim are in- 
troduced and the SEDs are likely to be better constrained when 
both data at 140 and 160 fim are added. When galaxies are only 
observed at 60, 100, and 140 /vm, a is found to be slightly larger 
(< Q'(measured) - Q'(predicted) >^ 0.2) than predicted from 
R(60,100) and the DH02 relation: we can tentatively explain this 
trend by the high values of a being likely to be ill-constrained 
without data at 160 jim leading to a degeneration of the proba- 
bility distribution and then to a slight overestimation of a. In all 
cases the difference in the a determinations remains lower than 
the uncertainty given by the Bayesian analysis. 

3.3.3. Constraints on the dust attenuation curve 

The presence of a bump at 220 run in the attenuation curve of 
external galaxies is still controversial. The best way to search 
for the presence of this bump is to work with UV spectra. In pi- 
oneering studies, Calzetti et al. ('1994 ) found no bump in local 



starburst galaxies, Burgarella et al. ( 2005 ) found some hints of a 
bump in one local Luminous IR Galaxy observed with GALEX, 
and, at higher z Noll et al. ( |2009a[ ) found evidence of a bump in 



star-forming galaxies with 1 < z < 2.5. Although the NUV band 
of GALEX lies on the bump itself ([Conroy etaD '2010), our SED 



fitting is unable to reliably ascertain whether a bump is present, 
as mentioned in section 3. We can say that the slope of the UV 
attenuation curve (parameter 6) is more reliable although the un- 
certainty remains large. The number of free parameters acting on 
the UV emission (star formation history, various parameters of 
dust attenuation) are at the origin of these poor determinations 



as discussed by 


Noll etal. (2009bl). 


Conroy et a 


. (2010|) adopted a more empirical approach and 



developed a diagnostic based on a colour analysis that clearly 
distinguishes models with or without a bump. The comparison 
of the FUV - NUV and NUV - u colours is sensitive to the 
presence of the bump. Their assumption of a constant star 
formation rate over 13.7 Gyr considerably reduces the number 
of free parameters. They also assumed a steep underlying 
curve (by adopting a low Ry) and found that their galaxies 
have FUV - NUV and NUV - u colours compatible with the 




Fig. 10. NUV - u versus FUV - NUV observed colours (black 
points) compared to Conroy et al. (2010 1 models (lines): the 
dashed blue line corresponds to a dust attenuation curve with 
no bump and the green solid line to a Milky Way bump strength. 



presence of a bump with an amplitude slightly lower than for the 
Milky Way extinction curve. In Fig [T O} we report the colours of 
our galaxies with the predictions of ,Conroy et al. (2010 1 for no 
bump or a bump with an amplitude similar to that of the Milky 
Way. Our galaxies span a larger range of colours than those 
explored by Conroy et al.|(|20101) but they clearly lie between the 
two predicted lines, suggesting that a bump is present that has a 
lower amplitude than for the MW. In all cases, the analysis of 
the full SED of these galaxies allowing for variations in the star 
formation history and dust attenuation scenarios cannot reach a 
reliable conclusion about the presence of a bump. 
The parameter 6 quantifies the variations in the slope of the 
attenuation law, the reference being the starburst redden- 
ing curve of [Calzetti et al. (2000i. For the measurement of 
< 6 >- -0.13 + 0.07, the slopes of the attenuation law are 
found to be slightly steeper than the Calzetti law but S is not 
well constrained with the available data (Fig|5]). This illustrates 
the difliculty in constraining the wavelength dependence of 
the attenuation curve with broad-band data only. As explained 
in Noll et al. ('2009b), S corresponds to a variation in Ry 
originally equal to 4.05 for the Calzetti et al. law: < 6 >- -0.13 
corresponds to Ry = 3.6. 



3.4. Age and mass of the stellar component 

The 4000 A break (hereafter D4000) provides information on the 
age of the stellar populations. It is small for young stellar popu- 
lations and large for old galaxies ( [Kaufl^mann et al.| [2003a). Th e 



D4000 definition adopted here is that of [Balogh et al.[ ( [i999[ ). 
It is expected to be insensitive to dust attenuation. Our mod- 
els have a resolution that is sufficient to calculate the D4000 
break (jMarastonj [2005 ), the calculation being based on the un- 
reddened spectra and thus an intrinsic, dust-free D4000 measure- 
ment. D4000 is quite sensitive to metallicity and its variation 
with the age of the stellar population is also strongly dependent 
on the star formation history (Maraston 2005 1 ). Its smooth and 



continuous increase when stars become older makes it a use- 
ful indicator of the age of the predominant stellar population. 
Moreover, its extensive use in analyses of SDSS data sets (e.g. 
Kauffmann et al. [2003b[ ) has provided references for the distri- 
bution of the D4000 in the nearby universe. In Fig[TT[ we report 
the D4000 values as a function of the stellar mass, both quan- 
tities being outputs of CIGALE. The value of D4000 (i.e. the 
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Fig, 11. Strength of the dust free 4000 A break, D4000, plotted 
against the stellar mass of the galaxies: both values are output 
parameters of the SED fitting code, the standard deviations given 
by the Bayesian analysis are plotted as error bars for 1/4 of the 
sample. 



Fig. 13. V-band attenuation. Ay, plotted against V-band attenua- 
tion of the young stellar population alone, Ay^ysp- Filled symbols 
are for galaxies with D4000 < 1.3 and empty ones for objects 
with D4000 > 1.3. Standard deviations are plotted as error bars 
for 1/4 of the sample. Crosses are over-plotted on symbols cor- 
responding to galaxies with /an < 0.2. 
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Fig. 12. Histogram of the specific star formation rate (SSFR) for 
the whole sample (heavy solid black line), and per mass bin (for 
the mass ranges 9.5 < log(Mstar) < 10 (thin solid blue line), 10 < 
log(Mstar < 10.5) (dotted green fine), and 10.5 < log(Mstar) < 1 1 
(heavy solid red line). 



age of the stellar populations) increases with the stellar mass, 
as also found by Kauffmann et al.| p003b) . Most objects have 
a low D4000 (D4000<1.5), even for massive galaxies. The sit- 
uation is quite different for the SDSS galaxies that exhibit a bi- 
modal D4000 distribution (roughly separated at D4000 = 1.5), 
the more massive galaxies (Mstar > 10" Mq) having D4000 
> 1.5. Our galaxies have D4000 values characteristic of galaxies 
with 10' < Mstai < lO'" Me in the SDSS. 
Another interesting parameter is the specific star formation rate 
SSFR - SFR/Mstai which is a measure of the present star for- 
mation activity relative to the past one and gives a rough indica- 
tion of the star formation history in galaxies. The histograms of 
SSFR for the whole sample and different mass bins are plotted 
in Fig [12] As expected from their D4000 and stellar mass distri- 
butions, galaxies exhibit a unimodal distribution of SSFR with 
a tail towards high SSFR, the average value for the whole sam- 
ple being 10"'" yr"' . This distribution is similar to that found for 
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Fig. 14. Distribution of the ratio of V-band attenuation over 
V-band attenuation of the young stellar population alone, 
Av/Av,ysp- The empty histogram corresponds to galaxies with 
D4000 > 1.3, the black filled histogram to galaxies with Z)4ooo < 
1.3 and the grey filled histogram to galaxies with /an < 0.2. 



galaxies with stellar mass between lO' and 10'° Mq in the SDSS 
( jBrinchmann et al.]|2004[ ). No obvious trend is found with stellar 
mass, in contrast to the findings of Brinchmann et al. (2004 1 in 
the SDSS. 

We can therefore conclude that our selection, which is mainly 
based on a detection at 140 yum with AKARI, is biased towards 
massive galaxies (< log(Mstar) >= 10.4 + 0.6[Mq]) that are quite 
active in star formation relative to the SDSS galaxy sample: the 
SSFR and D4000 found for our sample are representative of less 
massive galaxies when the galaxies are selected using visible 
data. 
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4. Dust attenuation 

4. 1 . Time-dependent dust attenuation 

Young stellar population s may be more affected by dust atten- 
uation than older ones dCalzetti et al.| |2000[ [Chariot & Falll 
1000} [Panuzzo et aE] [2007)i. CIG ALE allows us to separately 



vary dust attenuation for the two stellar populations considered 
for the star formation history by adjusting the parameter f^u, the 
reduction factor applied to the dust attenuation in the V band of 
the old population (expressed in magnitudes). Although f^a is 
poorly constrained, the global dust attenuation Ay (for both old 
and young stellar populations) and that found for the young stel- 
lar population alone Ay.ysp are clearly determined. We consider 
here D4000 as an estimate of the luminosity-weighted age of the 
stellar population and we divide the sample into two sub-samples 
corresponding to D4000 lower or larger than 1.3. In Fig 13 the 
total V-band attenuation Ay is plotted against the V band attenu- 
ation of the young stellar population only, Ay.ysp ■ The quantities 
are correlated and their relation depends on the stellar popula- 
tion age. The galaxies with D4000 > 1.3 for which is found 
to be lower than 0.2 exhibit a very low Ay: the old stellar popu- 
lation contributes substantially to the stellar emission and under- 
goes dust attenuation that is much lower than the one affecting 
the young component. This is indicative of a low value for the 
global attenuation in the V band relative to that applied to the 
young stellar population. 

Histograms of Av/Ay.ysp are plotted in Fig [14] As expected, the 
ratio is found to be lower for galaxies witnD4000 > 1.3 : < 
AvMv.ySP >= 0.43 ±0.17 (0.49 ±0.1 3 for sources with /att > 0.2) 
than for galaxies with D4000 < 1.3: < Av/Ay.ysp >= 0.70±0.14. 
[Calzetti et"aL] ( |2000| l found a factor of 0.44 between the color 
excess of the emission lines (related to the very young stellar 
populations) and the underlying continuum (older stellar popu- 
lations). Chariot & Fall (2000) recommend reducing by a factor 
three the effective optical depth for stars older than 10^ years. 
In the same way Panuzzo et al. ( 2007| l defined different age- 
dependent scale-heights for the stellar population, using again 
a typical timescale of ~ 10^ years to define young stars. In the 
present work, the young and old stellar populations are defined 
differently from these previous studies: the young stellar popula- 
tion does not represent stars younger than lO^years but the stellar 
population from a second episode of star formation and the pa- 
rameter D4000 measures a mean age of the global stellar popu- 
lations dominating the optical emission. Although the age of the 
young stellar population alone (f2) is ill-constrained (cf. section 
3), the mean value of the output parameter t2 is ~ 0.25 - 0.3 Gyr, 
i.e. older than lO^years. Moreover, we compare the attenuation 
of the global (old+young) stellar population Ay, which also de- 
pends on the star formation history, to that of the young stellar 
population alone, Ay^ysp. We cannot therefore directly compare 
with the difference of color excess E(B-V) found by Calzetti et 
|ar](r2000) between the nebular gas emission lines and the under- 
lying stellar continuum, which correspond to more extreme ages 
for the stellar populations (a few 10^ years for the emission lines 
and a few Gyr for the underlying continuum of the Balmer lines). 
The recipe of Chariot & Fall ( 2000), which consists of reducing 
by a factor three the effective absorption after 10^ years, is re- 
lated to our parameter f^a, which is defined as the ratio of A v for 
the old to young stellar populations. As discussed in section 3, 
/att is ill-constrained for individual cases. Nevertheless, /att = 1 
is chosen as the best-fit model when running CIGALE (x^ analy- 
sis) for only 4 % of the galaxies and the Bayesian analysis gives 
/att > 0.7 for only seven sources (cf. Fig. 7]). The mean value of 
/att for the whole sample is 0.39± 0. 15, be ow the mean of the in- 



put values (equal to 0.5, cf. Table|2]l, which, for an ill-constrained 
parameter, suggests a true mean lower than the return value of 
0.39. Hence we can conclude that /att < 0.5 as already found 
by [Noll et aL| ( |2009b| l for SINGS galaxies. An attenuation factor 
lower than 0.5 as well as a global dust attenuation (Ay) lower 
than that of the young stellar population alone ( Ay ysp) provide 
strong support for an age-dependent process for dust attenuation 
in galaxies. Given the uncertainties in the exact determination of 
/att and the different timescales adopted to define young stellar 
populations, our results are consistent with those of [Calzetti et 
al. (2000) and [Chariot & Fall| ( [2000) l and an age-dependent dust 
attenuation has to be assumed when performing SED fitting. 

4.2. Afuv as a function ofLi^/Lpuy 

The Lir/Lfuv raticj^is known to be a very robust estimator of 
Afuv as long as galaxies form stars efficiently. Several calibra- 
tions can be found, the main difference coming from the star 



formation history adopted (e.g. Buat et al. 2005 1 Cortese et al. 



2008 ). In quiescent systems, a substantial part of dust heating is 
due to photons emitted by old stars which do not emit primarily 
in UV. In this case, the relations between Lir/Lfuv and Afuv de- 
rived for an active star formation overestimate dust attenuation in 
the UV range. CIGALE gives the global attenuation in the FUV 
(combination of dust attenuation of young and old stellar popu- 
lations), which is plotted in Fig. 



15 



as a function of Lir/Lfuv- 
A very tight relation is found, as expected since our galaxies are 
found to constitute a homogeneous population of star-forming 
galaxies in terms of SSFR and D4000. Applying a polynomial 
regression on our data yields: 

Afuv = 0.483 -H 0.812 y -H 0.373 / H- 0.299 - 0.106 / 

, where y = log(LiR/LFuv)- 

We can compare these estimates to other Afuv - Ar/^'Fuv rela- 
tions from the literature, obtained with different scenaros for the 
dust attenuation and star formation rates. Gordon et al.[ ( [2000| l 
showed that the relation between Afuv and Lir/Lfuv is insensi- 
tive to dust attenuation configurations. Further studies identified 
a stronger sensitivity to the star formation history. ,Meurer et al. 
( 1999| l assumed a young stellar population only, for constant star 
formation rate over 1 to 100 Myr, to derive a relation between 
Afuv and IRX {IRX is defined as Lfir/Li6oo, the ratio of the far- 
IR luminosity between 4 and 120 fim to the luminosity at 1600 
A).|Cortese et al. ( 2008| explored more realistic models for nor- 
mal galaxies with exponentially decreasing SFRs and various r, 
Buat et al. (2005) obtained a mean relation for a large number of 
star formation scenarios, from short bursts to constant star for- 
mation rates. These different relations are reported in Fig. 15 All 



of them are obtained by assuming a single stellar/dust configu- 
ration, without any time dependence of the attenuation whereas 
a different attenuation for the old and young stellar population 
is introduced in CIGALE. The relation of Meurer et al.[ ( [1999| l 
systematically overestimates the dust attenuation we find in this 
work. Our results are more compatible with the relation of Buat 
let al.| ([2005 |l at least for galaxies with dust attenuations larger 
than 2 mag. The Cortese et al. (2008 ) relation with t - 6.2 Gyr 
provides a good fit for sources with intermediate attenuation (be- 
tween 1 and 2 mag). In all cases the variations in these estimates 
remain modest there being a difference of at most 0.3 mag be- 
tween our polynomial regression and the relations of Meurer et 



al. ( 1999) and Buat et al. (2005) that reaches 0.7 mag with the 
Cortese et al.| ( |2008| ) calibration for large dust attenuations. 

' Lpuv = y-Ly at 153 nm expressed in solar units as Lir 
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Fig. 15. Apuv plotted as a function of Lir/Lfuv, the result of Fig. 17. Apuv plotted against FUV - r. The sample is restricted 



the polynomial regression being plotted as a red thick solid line. 
Models from Meurer et al.| ( |T9991 i (black thin solid Une ), |Buat~et| 
[aTl (|2005 ) (blue dotted hne), and |Cortese et al.| ( |2008] l with t= 
6.2 Gyr (green dashed line) are overplotted. Standard errors in 
Apuv are plotted for 1/4 of the objects. 



to 265 galaxies with an error in FUV - r that is smaller than 0.2 
mag, the linear regression is plotted as a dashed line. The small 
dots correspond to FUV - NUV reported on the x-axis (similar 
to Fig. 16 1 and are overplotted here for comparison (see text for 
details). 




Fig. 16. Apuv plotted against FUV - NUV. The sample is re- 
stricted to 267 galaxies with an error in FUV - NUV lower than 
0.2 mag, and the linear regression is plotted as a dashed black 
line. Galaxies with Lir > 10" Lq are plotted as empty red cir- 
cles. Relations from previous studies are overplotted: Salim et aT] 
( |2007^ (green dot-dashed line), Seibert et al. (2005) (blue dotted 
line), and Meurer et al. ( 1999]l (black solid line). Standard errors 
in Apuv are plotted for 1/4 of the objects. 



4.3. Apuv and UV-optical colours 

Dust attenuation diagnostics based on UV-optical data alone can 
be very useful for applyung a correction for dust attenuation to 
the emission of galaxies not observed in far IR. |Meurer et al. 
( |199 9 ) found a relation between the slope of the rest-frame UV 
continuum /3 (defined as /i(erg cm^^s^'nm"') oc for A > 120 
nm) and dust attenuation traced by Li^/Lpuv for local starburst 
galaxies observed by lUE and IRAS. This local starburst relation 
is now widely used to estimate dust attenuation in UV-selected 
galaxies. In the local universe, the FUV - NUV colour is used 
as a proxy of /3, even if the NUV band overlaps the bump area of 
the Milky Way dust extinction curve (e.g. |Conroy et al.||2010[ )). 



Johnson et al. ( 2007| l investigated the Unk between Lir/Lpuv 
(and therefore Apuv) and various colours for galaxies detected 
by GALEX, SDSS, and Spitzer. Their sample represents both the 
red and blue sequence of the local universe and they found a de- 
pendence of Lir/Lpuv on D4000, i.e. on the age of the stellar 
populations. As described in section 3, our sample is strongly 
biased towards active star-forming galaxies, so we expect to find 
a smaller dispersion in the Apuv-colour relation than that found 
by [Johnson et al] (g007 ) . 

We compare the dust attenuation parameter Apuv obtained 
with our SED fitting method and two colours: the UV-only 
colour FUV - NUV and the longer base-line colour FUV - r. 
We first consider the FUV-NUV colour and restrict our analysis 
to sources with an error in FUV - NUV that is smaller than 0.2 
mag . The values are reported in Fig 16 and exhibit a fairly weak 



correlation (correlation coefficient R=0.59), a linear regression 
indicating that: 

Apuv = 2.39 (±0.20)x(FUV-NUV)+0.89 (+0.13) (cr = 0.83). 

Several authors proposed a relation between Apuv and FUV - 
NUV. [Meurer et al.| ( [1999[ ) obtained a linear relation between 
the slope of the UV continuum /3 and Afuv for starburst galax- 
ies, Seib ert et al.[ ( [2005 ) found that the starburst relation over- 
estimates Apuv by ~ 0.5 mag for a sample of galaxies not 
necessarily starbursting. These relations are reported in Fig 16 
wheren when necessary /3 is computed from the FUV - NUV 
colour (J3 = 2.23{FUV - NUV) - 2 ([Kong et al.[ |2004| ). We 



find that the relation proposed by Seibert et al. 



results more closely than that of Meurer et al 



pOOSFfits our 
199^ but, in 



most cases, it overestimates the dust attenuation for a given 
FUV - NUV colour Sah m et al.| (|2007| derived a shallower 
relation by fitting UV-optical data only (without IR data) for 
normal blue galaxies. Their relation is also reported in Fig [16] 
and appears consistent with our own regression. In all cases the 
dispersion in the data leads to a standard deviation around the re- 
gression line of 0.8 mag. The steepness of the relation between 
FUV-NUV and Apuv rnake the derivation of a dust attenuation 
from the FUV - NUV color quite uncertain. This weakness of 
the FUV - NUV colour as a dust attenuation estimator is also 
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underlined by |Conroy et al.| ( |2010| l. Using our regression for- 
mula leads to an uncertainty of 0.5 mag on Apuv for an error 
of only 0.2 mag on the FUV - NUV colour Luminous infrared 



galaxies (LIRGs) with Lir > 10" L© (Fig |16| are found to be 



closer to the starburst law than the bulk of the sample as akeady 
found by |Takeuchi et al.| ( |20T0l ) and |Howell etaL] ( pOTO) l, but the 
dispersion in their distribution remains large. 

The long baseline FUV - r is now considered, where we 
again restrict the analysis to sources with an error smaller than 
0.2 mag in this colour. A weak correlation (R=0.68) is found 
between Apuv and FUV - r. It is clear from Fig 17 that the 
larger range of values obtained for FUV - r compared to that for 
FUV-NUV leads to a more reliable estimate of dust attenuation 
for a similar uncertainty in the colours. A linear regression gives 

Afuv = 0.65 (±0.04) X {FUV - r) + 0.04 (±0.16) (cr = 0.76). 

The RMS dispersion cr is large reflecting the dispersion 
in the correlation. Nevertheless, the relation between Apuv 
and FUV - r is much flatter than that between Apuv and 
FUV - NUV implying an error of only 0.1 mag in Apuv for an 
error of 0.2 mag in FUV — r. [Johnson et al.| ( |2007 l) proposed 
several relations between Apyv, NUV - r, and D4000. For 
consistency, we prefer to compare FUV - r rather than NUV - r 
with Apuy. Nevertheless, we have checked that their relation 
corresponding to D4000 = 1 .25 (close to the mean D4000 found 
for our sample) is compatible with our values of Apuv and 
NUV - r. 

Both colours give only a crude estimate of Apuv with an 
RMS dispersion of ~ 0.8 mag. Since we have considered galax- 
ies actively forming stars the FUV - r colour is found to be a 
more reliable proxy than FUV-NUV of dust attenuation in UV: 
only one UV band is necessary and the relation with Apuv is less 
sensitive to the photometric uncertainty than when FUV-NUV 
is considered. In the same way, Johnson et al. ( 2007| l proposed 
to use the NUV - 3.6;um colour. 



5. SFR calibrations 

When both UV and IR fluxes are measured, a very eflicient way 
to estimate SFR is to add the contributions of both types of emis- 
sion to the total SFR (e.g. Iglesias-Paramo et al. 2006- jElbaz'etl 
al. 2007[ l. One must also account for dust heating by old stars, 
which is not directly related to the current star formation and 
given by 

SFRto, = (1 - /7)SFRiR + SFRpuv, 



where r] represents the fraction of IR emission due to dust heat- 
ing by old stars, SFRir being calculated by converting the to- 
tal IR luminosity Lir into SFR assuming that all the light from 
stars is absorbed by dust and SFRpuv is calculated from the ob- 
served UV luminosity Lpuy. The relation between the IR and 
FUV luminosities and the corresponding SFRs are derived from 
population synthesis models assuming a particular initial mass 
function and star formation history. In this work, we use the 
SFR calibrations from 'Buat et al.| ( [2008| ) for a constant SFR over 
10** years (the typical duration of the UV emission to reach a 
steady state) and a Kroupa initial mass function (K roupa||2001| l 
: log(SFRiR)M3yr-i = log(LiR)L3 - 9.97 and log(SFRpuv)Moyr-i = 
log(Lpuv)Lo -9.69. We prefer to combine SFRs rather than lumi- 
nosities to avoid bolometric corrections in UV, which depend on 
the wavelength range assumed for the UV (Bell et al. 2005 1 the 
FUV being taken to be monochromatic (at 153 nm) and directly 
calibrated in terms of the SFR. 
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Fig. 18. Fraction of dust heating by old stars, 77, plotted as a 
function of Lir, galaxies with Lir > 10" L© are plotted as empty 
red circles. Standard errors in both output parameters given by 
CIGALE are plotted as error bars for 1/4 of the sample. 



In Fig [18] we report the variation in 77 as a function of Lir for 
our sample, < r] >- 0.17 ± 0.10 for the whole sample. A small 
trend is found as a function of Lir with < rj >- 0.14 ± 0.08 
for Lir > 10" Lq (26 galaxies). Adopting 77 - 0.17 leads to a 
very good agreement between SFRtot and the SFR estimated by 
CIGALE (Fig. 



19 1. This agreement implies that the timescales 
over which both SFRs are calculated are similar As explained 
above, the calculation of SFRtot assumes a constant SFR over 
10** years. The SFR deduced from the code is mainly governed 
by the second stellar population, which is produced at a constant 
rate over t2 years (cf. Fig.|2]i. Even if t2 is poorly constrained (cf. 
section 3.2), its median value for the sample is 250 Myr, con- 
firming that the SFR measured by CIGALE is consistent with a 
constant SFR over ~ 10^ 



years. 



The value of 77 we find here is lower than that found by Hirashita 
|et al.|p003"] l (77 = 0.4). Befl (2003) estimated the mean contri- 
bution of old stellar populations to Lir to be 32% ± 16% for 
Lir < 10" Lq and 9% ± 5% for Lir > 10" Lq, these values and 
those of the present work are compatible within their uncertain- 
ties. 

da Cunha et al.| ( 2010| l estimated the fraction ^ of total IR lumi- 
nosity emitted by dust in the diff'use ISM, stars being assumed to 
migrate in this ambient ISM 10^ years after their birth. This pa- 
rameter is not directly comparable to 77 for several reasons: dust 
in the ISM is also heated by the young stars and a fraction of 
its emission is related to the star formation. We note that we de- 
fine 77 by assuming a constant star formation rate over 10^ years 
and stars in the ISM younger than 10^ years are also contribu- 
tors to the current SFR. We therefore expect ff, to be larger than 
77. Indeed, da Cunha et al. (2010 1 found median values of be- 
tween 0.5 and 0.7 with a dependence on the stellar mass and the 
specific star formation rate. We found a very weak positive cor- 
relation between 77 and Msi„r and 77 strongly correlates with the 
SSFR (Fig. 20 1 confirming the results of da Cunha et al. ( 2010) l. 
This result is unsurprising since in galaxies that are very active in 
star formation the dust is mainly heated by newly formed stars. 



6. Conclusions 

We have analysed the broad-band SEDs of 363 nearby galaxies 
from the UV to the far-IR by combining AKARI, SDSS, and 
GALEX surveys. The sample is close to representing one based 
on a far-IR selection. We fitted the SEDs with the code CIGALE 
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Fig. 19. Star formation rates obtained from CIGALE (x-axis) 
versus the total star formation rates SFRtot calculated with rj - 
0.17 (y-axis). 



discussed the reliabiUty of the FUV-NUV and Ff/V-r colours 
in estimating dust attenuation at FUV in the absence of IR data 
and find that both give only a very crude estimate of Apuv- 
The FUV - r colour is more reliable because the linear regres- 
sion Ap\jy-FUV - r is not as steep as that between Apuv and 
FUV - NUV. 

The SFR given by CIGALE is compared to the combination of 
SFRs estimated by adding the IR and the UV observed contribu- 
tions. Both estimates agree remarkably well when a contribution 
of dust heating by old stars of 17% is assumed. 
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Fig. 20. 77 plotted as a function of the specific star formation rate 
(SSFR), where galaxies with Lir > 10" are plotted as empty 
red circles. 



that performs a Bayesian analysis to deduce parameters related 
to the dust attenuation and star formation of each galaxy. An 
old (exponentially decreasing SFR) and a young (constant SFR) 
stellar population are combined and the attenuated stellar emis- 
sion is related to the dust emission on the basis of energy con- 
servation. We have shown that the detailed star formation his- 
tory is ill-constrained even when considering 11 photometric 
bands, whereas global properties such as the stellar mass, the 
total dust emission, the star formation rate, the 4000 A break, or 
the amount of dust attenuation are reliably estimated. 
Our galaxies appear to be both massive and active in star for- 
mation. Different dust attenuations are found for the young and 
old stellar populations leading to a global dust attenuation in the 
V band Ay that is lower than that assumed for the young stellar 
population only, by a factor of 0.4 for galaxies with a 4000 A 
break, D4000 > 1.3, and a factor of 0.7 for galaxies with D4000 
< 1.3 (D4000 being estimated from the SED fitting process). We 
have been unable to constrain the presence of a bump at 220 nm 
but have found that a dust attenuation law that is slightly steeper 



than that proposed by Calzetti et al. ( 2000 1 more closely repro- 
duces our data. 

The dust attenuation in the FUV band of GALEX is very well 
constrained by the Lir/Lfuv ratio. A polynomial calibration is 
provided and compared to other already published relations. We 



References 

Abazajian K. N., Adelman-McCarthy, J.K., Agueros, M.A. et al., 2009, ApJS, 
182, 543 

Balogh, M. L.; Morris, S. L.; Yee, H. K. C; Carlberg, R. G.; Ellingson, E. 1999, 

ApJ, 527, 54 
Bell, E. 2003, ApJ 586, 794 

Bell, E.F; Papovich, C, Wolf, C. et al. 2005, ApJ 625, 23 
Boissier S.,Gil de Paz, A., Boselh, A. et al. 2007, ApJS, 173, 524 
Boquien, M., Calzetti, D. , Kennicutt, R. et al., 2009, ApJ. 706, 553 
Brinchmann, J., Chariot, S.; White, S. D. M.2004, MNRAS, 351, 1151 
Buat, v., Iglesias-Paramo, J.; Seibert, M. et al. 2005, ApJ, 619, L51 
Buat, v., Boissier, S., Burgarella, D. et al. 2008, A&A, 483, 107 
Burgarella, D.; Buat, V.; Iglesias-Paramo, J. 2005, MNRAS. 360, 1413 
Calzetti, D., Kinney, A.L., Storchi-Bergmann, T.1994, ApJ 429, 582 
Calzetti, D., Armus, L., Bohlin, R. C. et al. 2000, ApJ, 533, 682 
CardeUi, J.A., Clayton, G.C., Mathis, J.S. 1989, ApJ, 345, 245 
Chariot, S., Fall, M. 2000, ApJ 539, 718 
Chary, R., Elbaz, D. 2001, ApJ. 556, 652 

Conroy, C, Schiminovich, D., Blanton, M.R. 2010, ApJ, 718, 184 
Cortese C, Boselli, A., Franzetti, P 2008, MNRAS, 386, 1157 
da Cunha, E., Chariot, S., Elbaz, D. 2008, MNRAS, 388, 1595 
da Cunha, E., Eminian, C, Chariot, S., Blaizot, J. 2010, MNRAS 403, 1894 
Daddi E., Alexander, D.L., Dickinson, M. et al. 2007, ApJ, 670, 156 
Dale, D., A., Helou, G., Contursi, A., Silbermann, N.A., Kolhatkar, S. 2001, ApJ, 
549,215 

Dale, D. A. , Helou, G. 2002, ApJ, 576, 159 

Dale, D. A., Gil de Paz, A., Gordon, K.D. et al. 2007, ApJ, 655, 863 
Elbaz, D., Daddi, E. , Le Borgne, D. et al. 2007, A&A, 468, 33 
Giovannoh, E., Buat, V., Noll, S., Burgarella, D., Magnelli, B. 2011, A&A, 525, 
150 

Gordon, K. D., Clayton, G.C., Witt, A. N., Misselt, K.A. 2000, ApJ 533, 236 
Iglesias-Paramo, J., Buat, V., Takeuchi, T. T., et al. 2006, ApJS, 164, 38 
Iglesias-Paramo, J., Buat, V., Hernandez-Fernandez, J., et al. 2007, ApJ, 670, 
279 

Hirashita, H., Buat, V., & Inoue, A. K. 2003, A&A, 410, 83 
Inoue, A.K., Buat, V., Burgarella, D., Panuzzo, P., Takeuchi, T.,T., Iglesias- 
Paramo, J. 2006, A& A, 370, 380 
Johnson, B. D., Schiminovich, D., Seibert, M. et al. 2007, ApJS, 173, 377 
Kauffmann, G., Heckman, T.M., White, S. et al. 2003, MNRAS, 341, 33 
Kauffmann, G.. Heckman. T.M., White, S. et al. 2003, MNRAS, 341, 54 
Kawada. M., Baba, H., Barthel, PD. et al., 2007, PASJ, 59, 389 
Kennicutt, R.C., Armus, L., Bendo, G. et al. 2003, PASP 115, 928 
Kong X., Chariot S., Brinchmann J., Fall S. M., 2004, MNRAS, 349, 769 
Kroupa P, 2001. MNRAS, 322, 231 

Howell J. H.,Armus, L., Mazzarella, J.M. et al. 2010, ApJ, 715, 572 
Maraston, C. 2005, MNRAS 362, 799 

Martin D. C, Fanson, J., Schiminovich, D. et al. 2005, ApJ, 619, LI 
Meurer, G. R., Heckman, T. M., & Calzetti, D. 1999, ApJ, 521, 64 
Noll, S., Pierini, D., Cimatti, A. et al. 2009a, A&A 499, 69 
Noll S., Burgarella D., Giovannoli E. et al. 2009b, A&A, 507, 1793 
Panuzzo P., Granato G. L., Buat V., Inoue A. K., Silva L., Iglesias-Paramo J., 
Bressan A., 2007, MNRAS, 375, 640 



14 



V. Buat et al.: AKARI-SDSS-GALEX galaxies 



Popescu, C.C., Misiriotis, A., Kylafis, N.D., Tuffs, R.J., Fischera, J. 2000, A&A, 

362, 138 

Popescu, C.C., Tuffs, R.J, Dopita, MA., Fischera, J.„ Kylafis, N.D., Madore, 

B.F. 2011, A&A, 527, 109 
Salim, S. Rich, M., Chariot, S. et al. 2007, ApJS, 173, 267 
Saunders, W., Sutherland, W.J., Maddox, S.J. et al. 2000, MNRAS 317, 55 
Seibert, M., Martin, D. C, Heckman, TM. et al. 2005, ApJ 619, L55 
Siebenmorgen, R., Kriigel, E. 2007, A&A 461, 445 
Silva, L, Granato, G.L., Bressan, A., Danese, L., 1998, ApJ 509, 103 
Schlegel, D.,J., Finkbeiner, D.R, Davis, M. 1998, ApJ, 500, 525 
Takeuchi T. T, Buat V., Heinis S. et al. 2010, A&A, 514, A4 
Tuffs, R. Popescu, C, Volk, H.J., Kylafis, N.D., Dopita, M.A. 2004, A&A 419, 

821 

Veron-Cetty , M.P, Veron, P 2006, A&A 455, 773 
Walcher, J., Lamareille, F, Vergani, D. et al. 2008, A&A 491, 713 
Walcher, J., Groves, B., Budavari, T, Dale,D. 2011, Ap&SS 331, 1 
Yamamura, I., Makiuti, S., Ikeda et al. 2010, AKARI-FIS Bright Source 
Catalogue Release note Version 1.0 



